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Abstract 

Gamma-ray burst statistics are best explained by a source population at cosmological 
distances, while spectroscopy and intensity histories of some individual bursts imply an 
origin on Galactic neutron stars. To resolve this inconsistency I suggest the presence of 
two populations, one at cosmological distances and the other Galactic. I build on ideas 
of Shemi and Piran (1990) and of Meszaros and Rees (1993) involving the interaction of 
fireball debris with surrounding clouds to explain the observed intensity histories in bursts 
at cosmological distances. The distances to the Galactic population are undetermined 
because they are too few to affect the statistics of intensity and direction; I explain them 
as resulting from magnetic reconnection in neutron star magnetospheres. 
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1. Introduction 

Attempts to explain all the observed gamma-ray bursts (GRB) with a single popula- 
tion of sources have become progressively more difficult. On one hand, their distribution 
on the sky has been observed, with steadily improving precision (Atteia, et al. 1987; Mee- 
gan, et al. 1992), to be isotropic, an observation which is naturally explained (Usov and 
Chibisov 1975; Goodman 1986, Paczyhski 1986, Mao and Paczyhski 1992, Fenimore et al. 
1992, Piran 1992) if they are at cosmological distances. On the other hand, a number of 
GRB have been reported to show (Higdon and Lingenfelter 1990) spectral features at a few 
tens of KeV and at about 400 KeV, which are readily interpretable as cyclotron lines and 
the two-photon positron annihilation line from the surface of magnetized neutron stars at 
Galactic distances, but which are inexplicable at cosmological distances. If their validity 
is accepted, the data appear irreconcilable. 

The problem is further complicated by the fact that straightforward models of radia- 
tion transport in GRB at cosmological distances (Goodman 1986; Paczyhski 1986) predict 
very brief bursts of radiation with thermalized spectra, in contradiction to observation, 
while attempts (Brainerd 1992; Katz 1992) to explain the spatial anisotropy and logA^ 
vs. log S or V/V max distributions of GRB in Galactic models require the assumption of a 
spherically symmetric halo of ~ 100 Kpc radius. Finally, the soft gamma repeaters (SGR) 
introduce additional confusion. The fact of their repetition and the identification of one of 
them (March 5, 1979; Cline 1980) with a supernova remnant in the LMC point strongly to 
a Galactic population, while the presence of spectral features and an 8 second periodicity 
(March 5, 1979) indicate origin on a magnetic neutron star. However, it is unclear that 
SGR should be considered GRB at all because their properties, including their spectra, are 
very different, and arguments made for SGR may be irrelevant to the problems of GRB. 

As a first step towards resolving the apparent inconsistencies, I consider the obvious 
possibility that there are two distinct populations of GRB. A cosmological population C 
accounts for most GRB, and explains the statistics of isotropy and the logA^ vs. log 5 or 
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V/V max distributions. A population G includes those GRB with spectral lines (a minor- 
ity), and originate on neutron stars at Galactic distances. An individual GRB cannot be 
assigned to a population unless it shows spectral lines, but the majority (probably over- 
whelming) of those without spectral lines must be members of population C. The SGR 
may be members of population G. I do not consider more exotic possibilities, such as GRB 
arising within the Oort cloud (Ruderman 1975), because they all seem far-fetched, even 
though I myself have discussed one of them (Katz 1993). 

In the absence of information about the intensity statistics and angular distribution of 
population G alone, it is not possible to discriminate between a disc and a halo origin. It 
will be important to obtain this information, which may be reducible from archival data. 

In §2 I discuss a possible mechanism for GRB at cosmological distances, building on 
recent suggestions by others. In §3 I more briefly discuss magnetic reconnection models of 
GRB at Galactic distances. The March 5, 1979 event in the LMC poses an acute problem of 
gamma-gamma pair production (Carrigan and Katz 1992) which must be faced, whether 
or not GRB of Population G have comparable distances and luminosities, and even if 
there is no population G. §4 contains a summary discussion. Unfortunately, unambiguous 
observational tests of the ideas discussed here will not be easy. 

2. Population C: GRB at Cosmological Distances 

The well-known failure of straightforward fireball models to explain the spectral and 
temporal properties of GRB led Shemi and Piran (1990) to consider neutrino-produced 
fireballs loaded with small (but not zero) amounts of ordinary matter; they found that 
the fireball could couple nearly all of its energy to the matter and (with the right values 
of the parameters) could accelerate it to relativistic velocity. Meszaros and Rees (1993) 
then pointed out that the interaction of this relativistic debris with surrounding matter 
might be characterized by times consistent with the range of GRB rise times and durations 
(10" 3 -10 3 sec). 
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I build on these ideas. The environments of GRB's at cosmological distances are open 
to much speculation (for example, are they low density galactic halos or dense nuclei of 
galaxies?), but the strong dumpiness of interstellar matter is a consequence of immutable 
atomic physics (cooling rates), and isolated discrete clouds are likely under a very wide 
range of conditions. The rarity of GRB makes it possible to assume favorable conditions, 
if these lie in the range of plausibility; there is no great difficulty if a considerably larger 
number of fireballs occurring in less favorable circumstances do not produce observable 
GRB. 

GRB at cosmological distances require the radiation of ~ 10 51 erg in observable 
gamma-rays. The complex chain of processes which lead to gamma-ray emission must 
be moderately efficient when the parameters do have favorable values, because the energy 
radiated as neutrinos by neutron star collapse, formation, or coalescence is unlikely to 
exceed 0.3M©c 2 ~ 5 x 10 53 erg, and may be considerably less. Thus, while we are enti- 
tled to assume favorable circumstances to explain the rare observable GRB, when these 
circumstances occur the resulting processes must be reasonably efficient. The fraction 
of neutrino energy converted to an electromagnetic fireball is small. Efficient conversion 
requires neutrino-neutrino collisions at angles in excess of 90°, but the neutrinos are gen- 
erally expanding outward in a neutrino fireball, with velocity vectors which are tending 
toward radial outflow. Optimal head-on collisions are particularly rare. The conversion 
of electromagnetic energy to particle kinetic energy also has an efficiency < 1. The final 
conversion to observable gamma-rays is the hardest part of the problem; it is easy to see 
how this could fail entirely. 

Relativistic invariants alone limit the amount of kinetic energy available for radiation 
by fireball debris. If a relativistic debris cloud with speed (3pc, Lorentz factor 7_p and 
proper mass per unit area a sweeps up a proper mass per unit area ctcx, then the efficiency 
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of radiation of the debris kinetic energy can be as large as 



a 



(1) 



e = 



a + 7 F (l - Pf)' 



Values of e > 1/2 are obtained for a > l/[^p(l + (3p)] m l/(2"f F ). Collisions with a very- 
broad range of clouds of circum-fireball matter are consistent with efficient conversion of 
kinetic energy to radiation. This is fortunate, because efficient production of GRB requires 
that common circum-fireball environments produce observable GRB in most directions; it 
is not possible to insist on fortuitous geometries or on special values of the parameters. 

Most of the kinetic energy will become available when the debris has swept up only 
a very little matter. If the energy of the explosion is Y then the proper mass of debris is 
Y/'Jfc 2 , and half the kinetic energy will become available when the swept- up proper mass 
is Y/2-fpC 2 . In a uniform medium of density p this will occur at an interaction radius 



where the numerical estimate assumed Y = 10 52 erg, p = 10~ 24 g/cm 3 and 7^ = 10 4 . 

The hardest part of the problem is turning the kinetic energy of the relativistic debris 
into the observed gamma-rays. The collision length of relativistic protons in ordinary 
matter is about 50 g/cm 2 , or about 100 Mpc at typical interstellar densities. Clearly, some 
collective process is necessary, and it must couple the proton and ion energy into that of 
electrons, which radiate more readily. Even relativistic electrons do not radiate rapidly 
under interstellar conditions; the radiation length of a 10 13 eV electron (corresponding 
to equipartition with a 7^ ~ 10 4 proton) for Compton scattering on a 3°K black body 
radiation field is ~ 10 23 cm, excessive by many orders of magnitude. 

In order to obtain short pulses of radiation at distances of order those given by (2) 
it is necessary that a coherent relativistically expanding front of radiating particles be 
directed nearly toward the observer. It is not sufficient that individual particles be observed 
only when directed towards the observer, a condition met by most relativistic radiation 
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processes. Therefore, ambient magnetic fields must not deflect particles significantly from 
their initial spherical expansion. This condition will be satisfied if the magnetic energy 
E mag in the interaction sphere of radius 77 is very much less than the debris energy, 
so the debris can sweep away the ambient magnetic field without significant deflection. 
If equipartition is assumed between the ambient magnetic field Bism and an ambient 
turbulent velocity field vism, then E mag /Y ~ v ]sm/1f c2 ^ V7f> so ^ ne ambient field 
may safely be ignored. 

a) Shock Structure 

When the debris shell collides with a cloud of ambient matter the resulting flow may 
be complex. If the shell and the cloud each initially had uniform density and velocity and 
negligible (on a relativistic scale) temperature, the geometry is slab-symmetric, and all 
bulk velocities are normal to the planes of symmetry, then the resulting shock structure is 
shown in Figure 1. There are two shocks SI and S2 and, in general, a contact discontinuity 
CD separating shocked fireball debris from shocked cloud. The equations relating the 
conditions in the four regions are cumbersome, except in the special symmetric case in 
which debris and cloud initially had the same composition and proper density. In this 
case, which I assume, there is no contact discontinuity and conditions in regions 2 and 3 
are identical, as are those in regions 1 and 4. 

The relativistic shock conditions (Landau and Lifshitz 1959) may be used to determine 
physical conditions. The thermodynamic variables are the proper internal energy density 
e, the proper pressure p, and the proper enthalpy density w = e+p. In the unshocked cloud 
ei = nim a c 2 , where n\ is the proper atomic number density and m a is the proper mass per 
atom, and p\ = 0. In the shocked cloud P2 = eil^ 3> e\ in the extreme- relativistic (ER) 
limit. I shall refer to the frame of the unshocked interstellar material as the local observer's 
frame; transformation to our frame requires application of the cosmological redshift. Then, 
to lowest nontrivial order in ei/e2 <C 1, the velocities of the fluids with respect to the shock 
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front SI are 



and" 1 



Vl 

c 



(gg ~Pi)(e2 

(e 2 - ei)(ei + p 2 ) 



1/2 



1-^ 

e 2 



W2 

c 



(P2 -Pi)(ei +P2) 



1/2 



3 V e 2 



(3) 



(4) 



(e 2 - ei)(e 2 +pi) 

The velocity discontinuity i> i 2 between fluids 1 and 2, measured in the frame of either, is 
obtained from the relativistic expression for the subtraction of velocities 



V12 ui/c-u 2 /c ei 
— 1 — 2 — 



1 — V1V2/C 2 



e2 



(5) 



this velocity is also the velocity u 2 l of shocked fluid 2 in the local observer's frame. The 
velocity vi is also the speed of the shock SI in that frame. 

Fluids 2 and 3 have the same velocity and, given our assumptions that n\ = 77.4 and 
ei = e 4 , the same values of the thermodynamic variables. Then the velocity of fluid 3 
in the frame of shock S2 is — v 2 , and the speed of fluid 4 in that same frame is —v±. 
The expressions for combinations of relativistic velocities may then be used to obtain the 
following results in the local observer's frame: 



^ = ^^1-2^ 

c c e 2 



c e 2 



1-2 



ei 



(6) 
(7) 
(8) 



<• \e 2 

It is now possible to calculate e 2 from the debris Lorentz factor 7^, defined in the 
local observer's frame, using Equation 8: 



If = 



e 2 



(1 - (Vc) 2 ) 1 / 2 2ei' 



(9) 



* Note the assertion in the first edition of Landau and Lifshitz that in the ER limit 
V2 — > c/3 1 / 2 is a typographical error; the correct limit w 2 — > c/3 is given in later editions. 
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e 2 « 2n 1 m a c 2 ^ F . (10) 

The Lorentz factor 72L of the shocked material in the local observer's frame is obtained 
from V2L = f 12, Equations 5 and 10: 

72L ~ If ■ ( n ) 

The detailed mechanics of the shock are obscure, but must be collisionless in order 
to form a shock at all. The shocked matter need not be in thermodynamic equilibrium. 
Heating of the shocked material by plasma instabilities is the source of dissipation; the 
distribution functions of particle energies will not be (relativistic) Maxwellians, but are 
more likely to be power laws. The distribution of energy between electrons and ions 
in uncertain. In a highly relativistic shock, as we expect here, electrons and ions are 
kinematically very similar (identical in the ER limit), so I will assume that the distribution 
of particle energies is independent of species; roughly half the post-shock energy resides 
in electrons. Any neutral matter does not interact with the shock, so that the density 
n\ refers only to the ionized component. At distances ~ 77 (Equation 2) the interstellar 
material will largely have been ionized by the flash of radiation associated with the fireball 
or by collision with debris. 

Using the shock jump conditions for the proper enthalpy w, wi = nim a c 2 , and the 
ER limit u; 2 ~ 4e2/3 yields the proper atomic density 

n2x2 (^X'\ (12) 
\m a c 2 J 



Define 72 by the relation 



72 = — (13) 



then the mean energy (in the frame of fluid 2) per particle is 72//m p c 2 , where \im v is the 
mean proper mass per particle. For pure ionized hydrogen \i = 0.5, while for the usual 
cosmic abundances (fully ionized) \i = 0.62. The mean Lorentz factor of an electron (in 



the frame of fluid 2) is 



72e=72 • (14) 

m P 



The proper density 77-2 (Equation 12) may be rewritten, using Equation 13, as 

n 2 « 4ni72, (15) 

reproducing the result 77-2 = 4ni for a strong but nonrelativistic shock (72 — > 1) in a gas 
with an adiabatic exponent of 5/3. The Lorentz factor 72 is found from its definition 
(Equation 13), and Equations 9 and 15: 



-»"2{£) ~{f) ■ (16) 

In the local observer's frame most of the particles are narrowly collimated in the direction of 
the motion of fluid 2, and the typical Lorentz factor is larger than those given in Equations 
14 and 16 by a factor ~ 72L (Equation 11). The angular width of collimation depends on 
the angular distribution of the particle momenta in the proper frame of fluid 2, which is 
unknown. If this is isotropic, then the locally observed angular width (for electrons as well 
as ions) 

do ~ 72L 1 « 7f 1/2 - (17) 

Note that this is a much broader angular distribution than the locally observed radia- 
tion pattern from a single particle, whose Lorentz factor is ~ (ions) or ~ ^Fm p /m e 
(electrons). 

b) Time Dependence 

The geometry of radiation from an advancing spherical shock is shown in Figure 2. 
The distant (but cosmologically local) observer may first see a flash of radiation from the 
fireball itself, whose arrival time is taken as t = 0. If the fireball is a consequence of 
the merger of binary neutron stars, as often assumed, the initial pulse includes bursts of 
neutrino and gravitational radiation, as well as electromagnetic radiation. Their emission is 
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essentially simultaneous, although their arrival may be affected by dispersion arising from 
plasma refraction, neutrino rest mass (if they have any), etc. The initial electromagnetic 
flash is expected to have a thermal spectrum and to be extremely brief (<C 10 -4 sec) 
because of the small size of the fireball (Goodman 1986; Carrigan and Katz 1992); if, as 
assumed here, the fireball energy is largely converted (Shemi and Piran 1990) to kinetic 
energy of debris this initial flash may be unobservably faint. However, if it is observed 
the time interval between it and the rest of the GRB is an important constraint on the 
emission geometry. 

Radiation emitted from a point (r, 9) on the expanding spherical shell arrives at the 
observer at a time 

fw r ( i-«»») +r /i_n 

\v cj (lg) 

r(l — cos 9) r 

~ 1 7 

c u 

where v is the shell's expansion velocity and the parameter (dimensionally but not physi- 
cally a velocity) u = vc/ (c — v). If the angular distribution of radiated intensity, measured 
in the local observer's frame, is /(#'), where 9' is the angle from the normal to the radiating 
surface, then the energy dE radiated by a patch of area dA is 

dE = f(9') dA. (19) 

Radiation directed toward the observer has 9' = 9. Using dA = 2nr 2 sin 9 d9 and dt = 
r sin 9 d9/c yields the observed power 

dE 

P(t,r) = — = 2ncrf(9). (20) 

The function f(9) is proportional to a convolution of the angular distribution of the ra- 
diating particles and their radiation pattern; as previously discussed, the latter is expected 
to be narrower than the former. A plausible guess is then 

fW - (21) 
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— 1/2 

where the angular width 9 ~ ^ F 1 is essentially the same as that of the momentum 
distribution (Equation 17). The pulse shape is then obtained from Equations 20 and 21, 
using Equation 18 to eliminate 9: 



P(t,r) oc 



t--)-+0 2 Q 

u/ r 



(22) 



0, t<£ 



where the approximation cos 9 l — 9 2 /2 has been used. P(£, r) has been plotted in Figure 
3, where the dimensionless parameter r = 2ctjr9\ has been defined. 

The pulse form of Equation 22 should be regarded only as an envelope, for the actual 
pulse shape will be modulated by the spatial distribution of matter which the debris shell 
sweeps up. One striking feature of Equation 22 is its abrupt rise, consistent with the 
observed rapid rises of some GRB. The characteristic width of this function is 

At~^~— . (23) 
2c 7 F c v ' 

Use of r ~ 77 (Equation 2) and ~ 10 4 yields At ~ 10 sec, the right order of magnitude 
for the duration of GRB. Much longer or shorter At may be possible for plausibly different 
values of the parameters, particularly the cloud density, which is uncertain even to order 
of magnitude. 

The debris shell and shock propagate into a very heterogeneous medium. The effects 
of structure in 9 are shown by the dashed lines in Figure 3, which assume a cloud uniform 
in the range 9\ < 9 < 92, with abrupt boundaries. A more realistic gradual density profile 
or shape would produce a gradual rise and decay; the abrupt rise remains if (and only if) 
the cloud includes the line 9 = 0. Thus this abrupt rise is expected for some, but perhaps 
not all, GRB, in accord with observations. 

A complete intensity profile of a GRB requires the integration of Equation (22) over 

r: 

put 

V(t)= / P(t,r)g(r)dr, (24) 
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where the weighting function g(r) includes both the fact that the energy available for 
radiation falls off for r > rj and the effects of dumpiness of the ambient matter as a 
function of r. The observable region is a narrow half-cone of apical angle ~ 6q: unless 
the scale of spatial structure is < 776*0 ~ 2 x 10 13 cm, dumpiness is more likely to be 
apparent as a function of r than of 9, justifying the use here of Equation 22 which ignored 
any dependence of density on 9. 

In the absence of spatial heterogeneity g(r) may be taken to impose a cutoff at r « 77, 
so that 

/•min(ri,ut) /.min(rj,ut) 2 J 

Pit) "i nw*J B ^TTFfW (25) 

The integral is elementary, but cumbersome, and of limited quantitative interest because 
of the artificiality of the assumed uniform density; the rise time is 77/11. Two possibilities 
should be distinguished: 

1. For a shock SI propagating through a homogeneous medium v/c ~ 1 — e\je^ ~ 
1 — l/(27i?) and u 2^pc Then, for 77 given by Equation 2 and 7_f = 10 4 , the rise 
time is several seconds, given by Equation 23 and inconsistent with very rapid rise 
times. This corresponds to a long GRB pulse envelope. 

2. On the other hand, the fireball debris propagates through a vacuum or very low density 
intercloud medium with a speed v/c = (1 — I/7I,) 1 / 2 , so that u m 2^pc; impacts upon 
small discrete clouds scattered within a region of size ~ 77 will only introduce a time- 
width ~ 10~ 3 sec, or less. This might not measurably broaden the abrupt rise given 
by Equation 22 and in Figure 3. 

Integration over r introduces two broadenings, one 0(ri/2^pc) associated with the 
entire emission region of size ~ 77, where the low intercloud density is appropriate, and 
another 0(r c /2^Fc) associated with individual clouds of size r c <C 77. These broadenings 
may each be much less than the envelope width (Equation 23), permitting an observed 
signal resembling that of Figure 3. Several sub-pulses may be observed if the debris shell 
collides with several isolated clouds in a medium dense enough to slow the intercloud shock, 
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so that interpulse times are 0(ri/2 r y F c). 

The actual situation is much more complicated than can be discussed here. For 
example, debris shell and clouds need not have the same proper densities, and each is 
likely to be spatially heterogeneous. Shocks propagating through heterogeneous media 
vary their strength, and produce associated continuous rarefactions and compressions. 
It is plausible that the complex structure of observed GRB could be explained by the 
interaction of relativistic debris shells with clumpy media, but more quantitative results 
would require numerical relativistic hydrodynamic calculations. 

c) Radiation 

The hardest part of the problem is turning electron kinetic energy into the observed 
radiation. Even though the typical Lorentz factor of an electron in the local observer's 
frame is ~ ^2Ll2e ~ A*7-F m p/ m e (from Equations 11, 14, and 16), their rate of synchrotron 
radiation and Compton scattering in plausible interstellar magnetic and radiation fields is 
low. I therefore make the radical suggestion that the collisionless shock produces approx- 
imate equipartition between the magnetic energy density and the particle energy density. 
The proper magnetic field B<i and energy density in fluid 2 are then, using Equation 10, 



where ( < 1/2 is a phenomenological parameter describing the approach to equipartition. 
The synchrotron energy loss time for an electron with Lorentz factor given by Equation 
14, assuming no correlation between the direction of the electron momentum and that of 
the magnetic field, is then, taking pure hydrogen composition 



The radiating volume is moving toward the observer with a bulk Lorentz factor 72L (Equa- 
tion 11), so that application of a Lorentz transform yields the local observer's measured 



B 2 

_1 = £ e2 ~ 2<>i7 F m a c 2 , 



(26) 




(27) 
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radiation time 

tabs = l2L(t r2 - V 2L t r 2/c ) 

= tvihiL (28) 

-C 1 (r^)~V LI xlO 7 sec. 
For a plausible interstellar cloud density n\ > 1 cm -3 and 7^ « 10 5 , t & s may be a 
millisecond or less. This justifies the assumption, made implicitly in the discussion of 
GRB rise times and pulse lengths, that shock-accelerated electrons radiate instantaneously; 
properly, the pulse profiles predicted by Equations 24 and 25 should be convolved with a 
broadening function which includes the radiation time, and which has a width t b s in the 
local observer's frame. 

The characteristic frequency of synchrotron radiation, measured in the frame of fluid 
2, is obtained from standard expressions using Equations 14 and 26. For pure hydrogen 
the result is 

^ 2 ~ Cl/2 (l^) ^ 3x1011 sec_1 ' W 
while Lorentz transformation to the local observer's frame, using Equation 11, yields 

v ohs ~ C 1/2 (r-^) 1/2 7l 3 x 10 11 sec" 1 . (30) 
V 1 cm 13 / 

MeV photons may be observed for £ = 1/2 if n\ ~ 1 cm -3 and 7f ~ 4 x 10 4 , for example. 

It has also been observed (Fishman 1993) that many GRB, or subpulses within them, 
show a progressive spectral softening with time. This is qualitatively explained using 
Figure 2 and Equation 18. If the radiation field is isotropic in the frame of fluid 2 (as will 
be the case if the particle distribution and magnetic field directions are isotropic) and has 
a characteristic photon energy, then in the local observer's frame the spectral hardness 
above this characteristic spectral peak will be a decreasing function of 9. Higher frequency 
photons are preferentially observed from smaller values of 9, which arrive earlier in the 
burst or sub-pulse, while lower frequency photons are observed over a wider range of 9 and 
hence over a longer time. A quantitative prediction for the spectral evolution with time 
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could be made by numerical integration of the synchrotron emission function, but would 
depend on (uncertain) assumptions made regarding the energy and angular distributions 
of the radiating electrons. 

3. Population G: Galactic GRB 

GRB which show spectral features, typically around a few tens of KeV and at 400 
KeV, have long been identified with Galactic magnetic neutron stars and are inexplicable 
at cosmological distances. Their distances cannot be determined from available data, 
and could be < 100 pc, ~ 100 Kpc, or anything in between. The familiar arguments 
concerning the mechanisms of GRB at Galactic distances center on two issues: the source 
of energy, and the physical conditions in the emitting region. The problems are harder, 
the greater the assumed distances. The observation of the March 5, 1979 event at a likely 
distance of 55 Kpc (Cline 1980) forces the consideration of distances of that order, and 
of correspondingly high luminosities, even though it is unclear whether it (a SGR) was a 
member of the Population G of GRB, or represented a distinct third class of objects. 

The central problem of distant and luminous gamma-ray sources is gamma-gamma 
pair production (Cavallo and Rees 1978; Schmidt 1978; Katz 1982; Epstein 1985; Carrigan 
and Katz 1992). This process does not permit the escape of a large luminosity of MeV 
gamma-rays from a small region unless they are collimated, and thus excludes many models 
of GRB or SGR at Galactic halo or cosmological distances. It is well known that this 
problem is avoided in a collimated relativistic outflow of radiating matter, a consideration 
which led to the popularity of fireball models, in which an opaque cloud of radiation and 
pair gas adiabatically expands and cools until its particles' velocity vectors are collimated 
outward. However, fireballs are conspicuously incapable of producing low redshift (400 
KeV) annihilation lines, line features at tens of KeV, or the observed long and complex 
time structure. The interaction of fireballs with their environment may solve the temporal 
problem, as discussed in §2, but offers no hope of solving the spectral problem. The case 
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for magnetic neutron stars for GRB or SGR with spectral lines remains as strong as the 
data. 

It is usually assumed that the radiating region of a GRB or SGR in a neutron star 
model is dominated by pair plasma, with n+ 3> rii, where n+ and rii are the positron 
and ion densities, respectively. This assumption is made, in analogy to fireballs, even for 
non-fireball models, perhaps because the threatened gamma-gamma pair production catas- 
trophe seems a likely source of dense pair plasma, and because the observed annihilation 
line requires a source of positrons. However, the assumption of pair dominance may not 
be justified in non- fireball models, such as are required to explain Population G GRB. 
If sufficient gamma-ray collimation is present to avoid a gamma-gamma pair production 
catastrophe, then the production of pairs may be negligibly small. When the observation 
of annihilation radiation provides empirical evidence for the production of some positrons, 
it should be remembered that an observably narrow annihilation line requires temperatures 
< 50 KeV, and may be produced by a comparatively small number of positrons precipi- 
tated onto the cool neutron star surface; a hot pair plasma does not produce a recognizable 
annihilation line. 

If a pair plasma is not an expanding fireball, it must be trapped on magnetic field 
lines (Carrigan and Katz 1992). Gravitation is unimportant for pairs, so they fill a magne- 
tosphere (presumably of a magnetic neutron star). However, they quickly (in a free-flight 
time) precipitate onto the stellar surface, where they annihilate, because they more rapidly 
radiate their transverse momentum by the cyclotron process; even if the radiation density 
is sufficient to maintain most leptons in excited Landau (magnetic) states (a condition 
satisfied under only the most extreme conditions), their interaction with the radiation 
field destroys their transverse adiabatic invariant. In this magnetosp here-filling geometry 
the emergent radiation, by whatever process, is not collimated, and gamma-gamma pair 
production imposes its usual limits on the emergent flux of MeV gamma-rays. 

It may be more satisfactory to consider an electron-ion plasma with only a small 
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admixture of positrons (sufficient to produce the observed annihilation line), so that n + <C 
rii. An electric field may accelerate the electrons, which radiate by bremsstrahlung or by 
the cyclotron process after elastic scattering on the ions raises them to excited Landau 
states. Because most of the leptons are negative, they form a broadly collimated beam; if 
they are relativistic the resulting radiation is similarly collimated and there is no gamma- 
gamma pair production catastrophe or limit (other than a Planck function at an effective 
temperature characterizing the electron distribution function) on the emergent intensity. 

In contrast, an electric field acting on a pair gas heats it but imparts no net momentum 
to the leptons; two counter streaming beams of gamma-rays readily produce pairs, rapidly 
achieving equilibrium with them and limiting the emergent intensity. A minority admixture 
of positrons in an electron-ion plasma produces only a proportionately small countercurrent 
of gamma-rays to those produced by the electrons. This countercurrent removes an equal 
current of electron-produced gamma-rays by pair production, but the remaining electron- 
produced gamma-rays form a collimated beam and escape comparatively freely, suffering 
little or no (depending on the degree of collimation) gamma-gamma pair production. This 
is a consequence of the net momentum imparted to the lepton-photon system by the electric 
field, in analogy to the momentum imparted to a sector of a fireball by adiabatic expansion. 

It is possible to make simple rough estimates of the parameters of an electrically heated 
ion-electron sheet plasma, which might be the source region of a GRB of Population G, 
following Katz (1993). Consider a sheet of thickness L, be composed of positive ions of 
charge Z and density (n e = Zrii), have transverse optical depth r and temperature 
T, and radiate power per unit area P. Define the dimensionless power per unit area 
p = Ph 3 / (rrigC 6 ), dimensionless thickness £ = Lm e c 2 /e 2 and temperature t = fc^T/ (m e c 2 ), 
where these quantities have been scaled to values characteristic of a relativistic electron 
(or pair) gas. The characteristic radiant intensity m^c 6 /h 3 = 4.3 x 10 35 erg/(cm 2 sec) and 
length e 2 /(m e c 2 ) = 2.8 x 10~ 13 cm (the classical electron radius). The optical depth is 

t ~ n e a L, (31) 
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where the characteristic cross-section a = e 4 /(mgC 4 ). At nonrelativistic energies the 
appropriate cross-section is the Thomson cross-section 87ro"o/3, but at semi-relativistic 
energies of interest do may be a fair approximation. The observation of a nonthermal 
spectrum implies that r cannot much exceed unity, but r ~ 1 and r<l are each possible. 

The large field of a magnetized neutron star has a number of effects. It enters the 
argument of the effective Coulomb logarithm in collisional processes, typically reducing 
it to In A In ( fee Tm e c/ \%eB)) (Katz 1982). Both bulk motion and current flow are 
restricted to be parallel to the field lines, justifying the assumption of thin sheet geometry 
and making the field distribution nearly force free (J x B = 0). Perhaps most important, 
it means that any electron energy resulting from motion perpendicular to the field is 
immediately radiated. Even in conditions characteristic of the March 5, 1979 event at 
55 Kpc the radiation density is far below Planckian (for t ~ 1), so that electrons may 
be assumed to be in their ground magnetic state until collisionally excited, and then to 
radiate as if in vacuum. The radiation rate per unit area may therefore be estimated using 
standard expressions for elastic scattering: 



This expression may be inverted, using the definitions of p, £, t, and a = e 2 /{Tic), to give 




(32) 



n 2 e e 4 LZ In A 
m e c t 1 / 2 




(33) 



the characteristic density m^c 3 /h 3 = 1.8 x 10 



cm 



. An alternative expression for n e is 



obtained from Equation 31: 



n e ~ 




(34) 



Equating the expressions (33) and (34) yields a result for the thickness: 



r 2 Z In A 

pt 1 / 2 ^ 



(35) 
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Very roughly, p ~ 10~ 4 for the March 5, 1979 event, if in the LMC, and p ~ 10~ 7 for 
a typical observed GRB if at 100 Kpc distance, corresponding to L ~ 1.4r 2 cm and 
L ~ 1.4 x 10 3 r 2 cm, respectively, where Z = 26, t = 0.1, and A = 10 were taken. The 
densities are correspondingly high. 

If the radiating sheet is driven by magnetic reconnection, as is plausible, then the 
radiated power may be related to the electrical work done: 

P ~ a el E 2 L, (36) 

where E is the electric field (properly, its component parallel to B) and a nonrelativistic 
expression (Spitzer 1962) is used for the (cgs) electrical conductivity: 

/2\ 3/2 (k B T) 3 / 2 
a ei ~ 2 - 



\*7 me /2 e 2 £ZlnA ^ 
t 3/2 m P c 3 



£ZlnA e 2 ' 

where the parameter £ > 1 is a correction factor which allows for the possibility of anoma- 
lous (plasma instability) resistivity when current densities and electron drift velocities 
are large and for the decrease in conductivity when electron velocities approach c. The 
characteristic conductivity m e c 3 /e 2 = 1.1 x 10 23 sec -1 . Equations (36) and (37) may be 
combined with the definitions of the various dimensionless parameters to give the electric 
field, current density, and electron drift velocity: 

„ fpa 3 £Z In A 2 \ 1/2 . . 

E~ ( - ^ 3/2 n e m e c 2 ) ; (38) 

«* = ~ W ~ W (40) 

The characteristic current density m 3 c 7 /e 5 = 6.5 x 10 38 esu/(cm 2 sec). The drift velocity 
is thus comparable to the electron thermal velocity unless £ S> 1; ion-acoustic instability 
is likely unless £ > rrii/{Zm e ). 
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GRB mechanisms based on magnetic reconnection, as in Solar flares, have been qual- 
itatively discussed for many years (Ruderman 1975). They may explain GRB energetics 
and phenomenology (Katz 1982). Magnetic reconnection by sheet currents provides a nat- 
ural explanation of the electrically heated sheets discussed in this section. If this model 
is assumed, then Maxwell's equations provide an additional relation among j, L, and the 
magnetic field B, permitting further constraints to be placed on the parameters. Because 

— » — » — » 

j is parallel to B, the direction of B rotates across a sheet current without changing its 
magnitude B. If the total angle of rotation across a uniform current sheet is 7r, then 

Defining the usual characteristic magnetic field B c = m 2 e c i /{eh) = 4.4 x 10 13 gauss and 
the dimensionless parameter b = B/{AB C ), Equation 41 may be rewritten 



rb^/ 2 Z In A 
pta 2 



e~ — . (42) 



Equating this expression to equation (35) yields 

r ~ —rpr- (43) 

For plausible t ~ 1 and b ~ 0.02, r ~ 10~ 4 ^ 1//2 ; either the emission region is very optically 
thin or the resistivity is dominated by plasma wave scattering, and is far in excess of its 
independent particle value. Either or both of these possibilities is acceptable. 

The electric field (Equation 38) may be evaluated, using Equation 34 for n e and 
Equation 35 for i, and defining the characteristic electric field E c = mlc 3 /(eh) = B c = 
4.4 x 10 13 cgs: 

^ (piZ In A m 4 c 5x 1/2 



V ^ 3/2 n 3 / (44) 

If the energy release is driven by magnetic reconnection it is proper to use Equation 43 for 



r, yielding 



E ~ ^E c . (45) 
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At Galactic halo distances the brightest GRB may have E sufficient to produce vacuum 
breakdown into a pair gas (Smith and Epstein 1993), but in less intense or closer GRB 
this will not occur, and the resistively heated ion-electron plasma discussed here may be 
sufficient. On a microscopic level, of course, the power is determined by the magnetic 
field strength and configuration, and by the mechanisms of plasma resistivity which drive 
reconnection. 

The observed non-thermal gamma-ray spectrum of GRB requires the presence of a 
nonthermal distribution of "runaway" electrons, consistent with the large Vd r and collective 
processes discussed above. It may be relevant that values of t ~ 1, consistent with the 
radiation of the bulk of the power of GRB at ~ MeV energies, correspond to a maximum in 
the conductivity and therefore, under conditions of a fixed potential drop, to a maximum 
in the power dissipation. 

If Z ^> 1, a e i increases approximately linearly with n + in the range Zrii < n + < Z 2 rii 
because of the increasing density of charge carriers without a corresponding decrease in 
their scattering length. This is in contrast to the usual near-independence of density of 
a ei. Pair production thus may provide a natural thermostat at t ~ 1, reducing the power 
dissipated in regimes at the pair production threshold by increasing the conductivity under 
conditions of constant current. It is evident, of course, that observable cyclotron and 
annihilation lines require require much lower values of t, and are plausibly produced by 
energy and positrons precipitated on the dense surface layers of the neutron star, cooled 
by black-body radiation. 

4. Discussion 

The fundamental problem of GRB phenomenology is the apparent inconsistency be- 
tween their spatial distribution, which points strongly toward a cosmological origin, and 
the spectral features observed in some GRB, apparently inconsistent with such an origin. 
In this paper I have tried to reconcile these apparently contradictory data by assuming 
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two disjoint populations of GRB. Because the argument for a cosmo logical Population C 
is statistical, while the argument for a Galactic Population G is based on the observation 
of spectral lines from only a minority of GRB (perhaps from only a minority of that Pop- 
ulation), it is difficult to assign an individual GRB to either population, unless it shows 
spectral lines or is identified with another astronomical object (the SGR of March 5, 1979 
is the only good extant example of such an identification). 

Fortunately, the models discussed in this paper predict another potential distinguish- 
ing characteristic. GRB in Population C produce their radiation by the synchrotron pro- 
cess of relativistic electrons. In the magnetic field is ordered the radiation will be linearly 
polarized. GRB in Population G produce their radiation by the cyclotron process of semi- 
relativistic electrons, Coulomb scattered into excited Landau states. This radiation is 
elliptically polarized, with a substantial circular component. In contrast, radiation pro- 
duced by annihilation in a pair gas, such as the initial burst from a fireball or the trapped 
pair gas discussed by Carrigan and Katz (1992), is unpolarized. 

The predicted characteristic frequency of gamma-ray emission in fireball debris im- 
pact models (Equation 30) is fairly sensitive to 7jt. The value of 7_f depends on physical 
conditions within the fireball (Shemi and Piran 1990), and a wide range of 7^ would be 
expected, with a corresponding range in spectra. In particular, smaller values of 7^ would 
lead to X-ray, ultraviolet, or visible bursts, which should be searched for. Bursts of lower 
frequency radiation should have longer durations and smoother time histories (Equation 
28), and perhaps also lower efficiencies, as accelerated electrons undergo adiabatic expan- 
sion before they radiate their energy. 

I thank P. C. Joss and I. A. Smith for discussions. 
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Figure Captions 



Figure 1: Flow geometry in frame of contact discontinuity CD. SI and S2 are shocks and 
numbers denote regions of fluid. 

Figure 2: Emission geometry. 

Figure 3: Radiation pulse emitted at radius r, plotted for 6q = 20c/ 'u. Dashed lines 
denote rise and fall for a cloud sharply circumscribed between 6>i = (ti6q — 2c/ u) 1 / 2 and 
#2 = ( r 2^o ~~ 2C/-U) 1 / 2 . Note that abrupt rise at r = 2c/(6qu) is obtained for a uniformly 
filled medium, independent of the angular distribution of radiation, and does not require 
a cloud bounded in angle. Units of the ordinate are arbitrary. 
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